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Abstract 

The spreading of liquid drops on surfaces corrugated with micron-scale parallel grooves 
is studied both experimentally and numerically. Because of the surface patterning, the 
typical final drop shape is no longer spherical. The elongation direction can be either 
parallel or perpendicular to the direction of the grooves, depending on the initial drop 
conditions. We interpret this result as a consequence of both the anisotropy of the contact 
line movement over the surface and the difference in the motion of the advancing and 
receding contact lines. Parallel to the grooves, we find little hysteresis due to the surface 
patterning and that the average contact angle approximately conforms to Wenzel's law 
as long as the drop radius is much larger than the typical length scale of the grooves. 
Perpendicular to the grooves, the contact line can be pinned at the edges of the ridges 
leading to large contact angle hysteresis. 

1 Introduction 

In recent years wetting and spreading phenomena have received continued attention from 
the scientific community due to their broad application in, for example, microfluidic devices, 
surface coating and biomimetics. Surface roughness can be exploited to significantly alter 
the behaviour of fluids moving over a surface. Examples include plants where micron-scale 
bumps on the leaves lead to superhydrophobic behaviour [T], desert beetles who use hy- 
drophilic patches on their back to collect dew [2], and butterfly wings which are patterned 
anisotropically to promote directional run-off [3]. 

It is now possible to reproduce heterogeneous surface patterning in a very controlled manner 
on micron length scales. Regular arrays of chemical patches [H El E] and posts [3 El EJ [TO] 



1 



of different shapes and sizes are regularly fabricated and several authors [T2t [T3] have 
even shown the possibilities of manufacturing multi-scale surface patterns. Recently such 
patterning has been used to control the movement of drops |16[ [T7] and to attempt to enhance 
flow in microchannels \18\ [T9] . 

The motion of drops on patterned surfaces is complex because of pinning and hysteresis. In 
particular surfaces with anisotropic patterning on a scale comparable to the drop size can 
result in elongated drop shapes and in a different motion parallel and perpendicular to the 
grooves. Gleiche et al. [1] showed that there is an anisotropy in the average value of contact 
angle and contact angle hysteresis on chemically nanostructured surfaces. Brandon et al. |20j 
investigated the effect of drop size on chemically striped surfaces and found that the drop 
anisotropy and contact angle hysteresis depended on the drop volume. Elongated drop shapes 
were also obtained by Chen et al. [21] and Chung et al. [22] for hydrophobic and hydrophilic 
grooved surfaces respectively. Narhe and Beysens [23] studied the growth dynamics of water 
drops condensing on grooved surfaces and showed that similar elongated drop shapes can 
be found during growth when the surface is hydrophilic, but is absent when the surface is 
superhydrophobic. Pakkanen and Hirvi [24j further showed that the anisotropy still persists 
when surface patterning is nanoscopic. Morita et al. [5] and Yoshimitsu et al. [10] studied 
the dynamics of drops sliding on chemically striped and hydrophobic grooved surfaces and 
found that the sliding angles are considerably larger for drops moving perpendicular to the 
stripes. 

In this paper we present a number of experiments and lattice Boltzmann simulations of drops 
spreading or dewetting on a hydrophilic surface patterned with parallel grooves. Our results 
highlight the importance of hysteresis and energetic barriers due to the surface patterning: 
because of the asymmetric behaviour of the advancing and receding contact line, the final 
drop shape can be elongated either along or perpendicular to the grooves. When the contact 
line is advancing, the drop is found to be elongated parallel to the grooves. On the other 
hand, when the contact line is receding, the drop is elongated perpendicular to the grooves. 
We present both equilibrium and quasi-static drop experiments and simulations in section 
4. In section 5, an investigation of the influence of three dimensionless surface parameters is 
presented: the roughness factor, the aspect ratio of the barriers and the dimensions of the 
drop versus the dimensions of the barriers. As expected, the contact angle, particularly that 
of the interface lying perpendicular to the grooves, strongly depends on both the parameters 
of the surface and the history of the drop motion. The complementary nature of experiments 
and simulations is also demonstrated. Sections 2 and 3 summarise the experimental and 
numerical approaches respectively. 
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Figure 1: Periodically corrugated surfaces used in the (a) experiments and (b) simulations. 

2 Experiments 

Two sets of experiments were performed with anisotropically grooved surfaces. 

In the first set of experiments, micron-sized droplets were placed on sub micron-scale cor- 
rugated polyimide surfaces. The polyimide layers were created on glass by spincoating a 
commercially available mixture (JSR AL-3046) of polyimide from solvent. After spincoating, 
the samples were heated for 5 minutes at 80°C on a hotplate to evaporate the solvent, and 
afterwards thermally cured in a vacuum oven at 180°C for 1.5 hours. Afterwards the samples 
were rubbed unidirectionally with a velvet cloth with a force ~ 1 kg/cm^ over a distance of 
30 cm with a velocity of 10 cm/s. This created grooves that were determined by AFM to 
be between 15 to 20 nm deep and 50 to 200 nm wide. Non-rubbed surfaces were prepared 
identically, only the rubbing step was omitted. 

Surface energy measurements were made on the prepared polyimide surfaces. By measuring 
the contact angle of water, ethylene glycol and diiodomethane drops on the surface, the surface 
energy was calculated to be 43.6 mN/m, employing the OWRK method. 

Drops consisting of a monomer of 99.0 w% ethoxylated bisphenol-A dimethacrylate (Sartomer 
Europe BV), having a surface tension of 41 mN/m, together with 1.0 w% photo-initiator (Ir- 
gacure 184, Ciba Specialty Chemicals), were placed on the surface by two methods. In the 
first method, the mixture was sprayed manually onto the polyimide surface. The resulting 
drops were placed randomly on the surface and were of varying size. In the second method, 
micro-transfer printing [H], a PDMS elastomeric stamp with an array of trapezoidal posts 
with face areas of 100 x 100 /im was pressed into a thin spincoated layer of monomer mixture 
and afterwards pressed against the polyimide layer. Part of the mixture was transferred to 
the polyimide surface, after which the stamp was detached from the surface leaving an array 
of drops with equal size. The drops were polymerized within one minute after placement. 
This fixed the geometry of the drop and enabled analysis afterwards, while still allowing the 
drop sufficient time to equilibrate. The drop x — y — z profiles were determined by optical in- 
terferometry (Fogale Zoomsurf 3D), and the parallel and perpendicular cross-sectional height 



3 



profiles were extracted and ellipse-fitted to determine the contact angles and elongation. A 
minimum of five drops were examined for each surface. 

In the second set of experiments, smooth corrugated polymer surfaces were created by means 
of a four-step process called photo-embossing [T^. First a so-called photo-polymer layer 
was spin-coated on glass forming layers approximately 20 fim thick. The mixture consisted 
of 41.9 w% monomer dipentaerythritol penta/hexa-acrylate (Sigma Aldrich), 41.9 w% poly- 
mer polybenzylmethacrylate (Mw 70 kg/mol; Scientific Polymer Products), 4.20 w% photo- 
initiator Irgacure 819 (Ciba Specialty Chemicals) and 12.0 w% inhibitor tert-butyl hydro- 
quinone (Aldrich). In the second step, the layer was illuminated with UV light through a 
suitable lithographic mask (described below), partially polymerizing the exposed parts of 
the layer. In the third step, the sample was heated and part of the monomers in the non- 
illuminated regions diffused towards the illuminated regions, creating the corrugated surface. 
In the last step the surface was flood-exposed with UV light to polymerize the full layer. 
Similarly to the polyimide surface, the photo-embossed surface was found to have a surface 
energy of 44.7 mN/m. 

The lithographic masks employed had sizes of 1 x 1 cm and consisted of parallel stripes, where 
each stripe was as wide as the distance between the stripes, i.e. 50 % of the incident light is 
transmitted. The pitch, i.e. the repeat length of the mask (line to line) was varied from 10 fj,m 
to 100 fim. The depth of the grooves could be varied between 100 nm and 5 fj,m. by varying 
the illumination time, with the exception that the maximum depth could only be reached 
with masks with pitches above 40 /im because sufficient material needs to available to be be 
displaced by the photo-embossing process. The depth of the grooves was again determined 
with optical interferometry. 

Water drops of 4 fil were placed on the photo-embossed polymer surfaces by an automated 
dispensing needle (Optical Contact Angle setup, Dataphysics) at a constant height above the 
surface to minimize any variations in drop placement. The water used for the measurements 
was purified (Millipore Super-Q) and CO2 was removed (Elix-10 UV). The drop properties 
were measured on the surface with the contact angle measurement setup . Both and 6± 
were determined, as well as the drop elongation e. The contact angles and base radii were 
determined by means of digital drop shape fitting employing an ellipse fitting algorithm. 

The two sets of experiments are complementary, in that they employ different surfaces and 
probe liquids. The rubbed surfaces have nano-meter scale corrugations, whereas the photo- 
embossed surfaces have micron-scaled corrugations. Micro-transfer printing is suitable for 
creating well-defined dewetting drops, while photo-embossed surfaces have the great benefit 
of being smooth as well as controllable with regards to the size of the corrugations. In the 
following sections the results from both sets of experiments are presented next to each other 
highlighting different aspects of anisotropic hysteresis. 
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3 Lattice Boltzmann simulations 



We use a mesoscale simulation approach where the equihbrium properties of the drop are 
modelled by a continuum free energy 



'1' = J^iMn) + ^{dany)dV + j^^l:s{ns)dS. (1) 

ipbiiT-) is a bulk free energy term which we take to be [26] 

Mn) = Pcii^n + ifii^l - 2i/„ + 3 - 2/?r^) , (2) 

where = (n — nc)/nc, Tyj = (Tc — T) /T^ and n, ric, T, and pc are the local density, critical 
density, local temperature, critical temperature and critical pressure of the fluid respectively. 
This choice of free energy leads to two coexisting bulk phases (liquid and gas) of density 



nc(l =b \J I3tw)- The second term in Eq. ([T]) models the free energy associated with any 
interfaces in the system, n is related to the surface tension via 7 = {4:y/2Kpc{PTw)^^'^nc)/3 
[25] . The last term in Eq. ([T|) describes the interactions between the fluid and the solid 
surface. Following Cahn [27] the surface energy density is taken to be ips{n) = —(pUg, where 
Us is the value of the fluid density at the surface. The strength of interaction, and hence the 
local equilibrium contact angle, 9^, is parameterized by the variable (p. Here the simulation 
parameters are chosen to give nuquid = 4.11, ngas = 2.89, 7 = 5.14 x 10~^, and 9e = 70°. 
The typical drop and pattern sizes are of order 100 and 10 lattice spacings respectively. 
Simulation and physical parameters are related by choosing a length scale Iq, a time scale to 5 
and a mass scale tuq. A simulation parameter with dimensions [^]"^[i]"'^[m]"^ is multiplied 
by [^o]"'^[^o]"'^['TT'o]"^ to give the physical value [5H]. More specific details on the model used 
here, including the way we have implemented complicated surface geometries, can be found 
in [29]. 



The dynamics of the drop is described by the continuity ([3]) and the Navier-Stokes equations 

m 

dtn + da{nua) = , (3) 
dt{nua) + dpinuaup) = -dpPap + vd/ilnidpUa + d^Ufs + da^d^u^)] , (4) 

where u, P, and are the local velocity, pressure tensor, and kinematic viscosity respectively. 
The thermodynamic properties of the drop are input via the pressure tensor P which is 
calculated from the free energy [26^ [29l [28] via 

-Pq/3 = ipb{n) - ^{dauf - Knd^^n)6ai3 + K{do,n){di3n) , (5) 
Pb{n) = Pcii^n + l)H3^n -2un + l- 2/3r^). (6) 
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Eqs. §i) and (g]) are solved using a lattice Boltzmann algorithm [26 1 [29l[28l[30l [3T] . 

This model (or very similar approaches e.g. [32l [33l [S] ) has been shown to be a useful tool to 
study several aspects of drop dynamics on patterned surfaces, such as the problem of contact 
angle hysteresis [35] , and drops spreading on super hydrophobic [29] and chemically patterned 
surfaces [28] , 

To render the simulation feasible, the typical surface pattern used in the simulations is slightly 
different to that produced experimentally. It is therefore important to note here that our aim 
in this paper is to obtain a consistent qualitative understanding of the problem rather than 
to attempt to exactly match experiments and simulations. The shape of the grooves is taken 
to be rectangular, as shown in Fig. W[h). We define the aspect ratio as the ratio of the groove 
height to the pitch. No-slip boundary conditions are imposed on the velocity field on the 
surfaces adjacent to and opposite to the drop and periodic boundary conditions are used in 
the other two directions. 

4 Anisotropic drop morphology on a corrugated surface 

In this section, the anisotropy in the final shape of a drop spreading on a corrugated sur- 
face is investigated. First, we performed experiments where liquid drops are sprayed or 
micro-transfer printed onto the rubbed polyimide surface. We find that even though both 
preparation methods result in an elongated drop shape, the direction of elongation relative 
to the surface grooves is different. The drop is elongated parallel to the grooves if sprayed, 
while micro-transfer printed drops are elongated perpendicular to the grooves. 

Results consistent with the first set of experiments above were found when water is quasi- 
statically added to a spreading drop prepared on a photo-embossed surface. Here different 
qualitative behaviours are observed parallel and perpendicular to the grooves, further high- 
lighting the asymmetry of contact line motion in these two directions. The experiments are 
then complemented by performing lattice Boltzmann simulations to further understand the 
role of the surface anisotropy in determining the final shape of the drop. 

4.1 Experiments 

In Fig. Wljo) we show a typical final shape of a liquid drop that is sprayed onto the surface, 
causing the the drop contact line to advance to wet the corrugated surface. Since the sur- 
face patterning is not isotropic, the advancing contact line behaves differently parallel and 
perpendicular to the grooves. Based on Johnson and Dettre's work [36] no hysteresis due to 
the surface patterning is expected in the parallel direction as there are no energetic barriers 
present which could tend to pin the contact line. On the other hand, perpendicular to the 
grooves, surface undulations are known to pin the contact line [371 [38l I39j . Contact line 
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Figure 2: Typical drop shape observed in experiments when (a) there is no surface patterning, (b) the drop 
is sprayed onto the surface and spreads, and (c) the drop is stamped onto the surface and dewets. Although 
not visible in this image, the surfaces of (b) and (c) were patterned with sub micron-scale ridges. Here the 
ridges are oriented in the horizontal direction, (d), (e) and (f) Lattice Boltzmann simulation results when 
(d) the drop spreads on a smooth surface, (e) and (f) the drop volume is quasi-statically (e) increased and 
(f) decreased. The experimental drop elongation parallel (b, e) and perpendicular (c, f) to the grooves is 
clearly reproduced in the simulations, (g) Experimental measurements of the average elongation as function 
of printing temperature for the printed aery late drops on rubbed and non-rubbed polyimide. The drops on 
non-rubbed surfaces have approximately spherical contact lines at all temperatures, whereas with decreasing 
temperatures the drops on rubbed surfaces become more and more elongated. 
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Figure 3: Graphical illustration of Gibbs criteria: (a) advancing and (b) receding contact line. 

pinning leads to strong contact angle hysteresis and this causes the advancing contact angle, 
0A, to differ from that in the parallel direction. 

The pinning is illustrated in Fig. EI|a). For the drop contact line to advance, it has to wet the 
side of the grooves, which, according to the Gibb's criterion [38l|39], occurs when 6a = 9e+90° 
for rectangular ridges. More generally 9a = Oe + cx in two dimensions where a is the maximum 
inclination of the surface. In three dimensions, the value of the advancing (and receding - 
see below) angle is not extreme as that predicted by the Gibbs' criteria, due to the energy 
costs associated with the surface deformation from the spherical cap shape. The advancing 
angle is, nonetheless, generally larger than the advancing angle parallel to the grooves. As a 
consequence, it is easier for the drop to spread in the parallel direction and hence the drop 
shape is elongated parallel to the grooves as shown in Fig. [H[b). This is the final drop shape 
typically found in corrugated surface experiments |2H [2^ I22j . The shape of a drop on a 
smooth surface is also shown in Fig. [2]^a) for comparison. 

When the drop is (micro-transfer) printed onto the surface, on the other hand, the drop 
contact line retreats to dewet the surface, as the liquid is initially spread out much further 
than its equilibrium shape. This occurs when the drop contact angle is smaller than or equal 
to the receding contact angle. If no contact angle hysteresis is present, the final drop shape 
will be independent of the initial conditions. However, as discussed above, hysteresis is an 
important effect in the perpendicular direction and as a result, the receding angles are again 
different parallel and perpendicular to the grooves. Contact line pinning for the receding 
motion is illustrated in Fig. EJb). For the contact line to recede, the drop has to dewet 
the side of the posts, which happens when Or = 9(. — 90° for rectangular ridges. For other 
geometries, Oii = 9(. — am. two dimensions where a is the maximum inclination of the surface. 
Since the receding angle in the perpendicular direction is smaller than the receding angle 
parallel to the grooves, it is easier for the drop to dewet in the parallel direction and the final 
drop shape is typically elongated perpendicular to the grooves, as shown in Fig. [2] (c). 

The effect of the printing temperature on the final shape of the printed drops showed an 
interesting result, as shown in Fig. El^g). With increasing temperature, the elongation per- 
pendicular to the barriers decreases. Clearly, the pinning ability of the barriers is decreased 
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Figure 4: Advancing contact angle measurement for a drop of water spreading on a photo-embossed surface, 
(a) and (c) are measurements of the base diameter (in arbitrary units because the setup was not length- 
caUbrated), while (b) and (d) are the contact angles as a function of time, (a) and (b) are the diameter and 
contact angle perpendicular to the grooves, (c) and (d) show the same variables parallel to the grooves. The 
initial drop volume was 4 jA. Water was added with 0.06 /il per second effecting a quasi-static, linear increase 
in volume. The surface had a pitch 80 ^im and groove depth of 4.3 fim, i.e. an aspect ratio of 0.10. The 
equilibrium contact angle of the surface was 6e ~ 70°. The characteristic effect of crossing surface barriers in 
the perpendicular direction is clearly visible in (a) and (b), while (c) and (d) demonstrate that in the parallel 
direction the motion is much smoother. 
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at higher temperatures and the drop morphology is less affected by hysteresis. A possible 
explanation for this effect is that the viscosity at higher temperature decreases so that the 
drop in general and the contact line in particular move faster and dissipate less energy. Even 
without determining the exact barrier crossing mechanism, in general the interface can more 
easily overcome pinning on the ridges by having more energy available for deforming its sur- 
face locally, enabling it to cross more barriers. This is further supported by our simulation 
results, where we found similar dependence of drop elongation on the drop viscosity and 
impact velocity. The typical drop elongation decreased as we reduced the drop viscosity or 
increased the drop initial kinetic energy. 

The mechanism by which the interface crosses the ridges can be envisaged as the nucleation 
mechanism formulated by DeGennes [25]: in this process an advancing drop will first form 
a nucleus in the next grooves, after which the interface spreads along the grooves until a 
new and stable, but less elongated morphology is reached. For drops retracting over the 
barrier, the inverse process occurs. The drop retreats parallel to the grooves until it becomes 
sufficiently distorted that it is favourable for it to retract, dewetting the outermost grooves. 
We were able to observe this retracting mechanism with evaporating water drops which were 
placed on a photo-embossed surface. This is shown in Fig. [5l 

Our interpretation of the role of contact line pinning in determining the drop shape is backed 
up by the experimental evidence obtained when considering a quasi-statically growing drop 
on a photo-embossed surface, shown in Figs. HI The drop contact angles and base radii both 
parallel and perpendicular to the grooves are measured as a function of drop volume. Three 
distinct features are visible: (i) the saw-tooth shaped variation with time of the perpendic- 
ular contact angle, (ii) the step-shaped variation of the perpendicular base radius and (iii) 
the absence of these distinct features in the parallel direction where the base radius increases 
continuously and the contact angle is roughly constant over the whole range of the experi- 
ment. Features (i) and (ii) are clear indications of contact line pinning in the perpendicular 
direction, while feature (iii) shows that pinning of the contact line due to the surface pattern- 
ing does not occur in the parallel direction, leading to a significantly lower advancing contact 
angle. These features were reproduced for all aspect ratios and surfaces. For receding contact 
line, no depinning was observed in the parallel direction, while depinning in the perpendicular 
direction occurred irregularly, with interfaces sometimes crossing multiple ridges at once. For 
this reason, reproducibly quantifying the depinning as function of the surface parameters was 
unsuccessful. We found that, on the unpatterned surfaces, receding contact angle measure- 
ments give values of about 15° to 25° with occasional outliers of up to 50°. On the patterned 
surfaces, receding contact angle measurements in the parallel direction give similar values to 
those for the unpatterned samples, while the receding angles in the perpendicular direction 
are on average about 5° lower than those in the parallel direction. The high experimental 
scatter is attributed to the common problem of measuring very low receding contact angles. 
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Figure 5: (color online) Microscope observation of the shape of the receding contact line as a water drop 
evaporates over time. The drop retreats parallel to the grooves until it becomes sufficiently distorted that it 
is favourable for it to retract, dewetting the outermost grooves. 1% Rhodamine B is added to the drop for 
visibility. 
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Figure 6: Simulations of drop contact angles on a grooved surface as a function of the drop volume for 
groove height h — 3, ridge width wi = 10 and pitch p — 40. The drop was initialized at V ~ 0.75 x lO'' and 
equilibrated. Then the drop volume was either slowly decreased or increased. 

namely that incidental pinning and depinning are observed even on (apparently) smooth sur- 
faces. The base radii were also greater in the perpendicular direction, further proving that 
hysteresis is increased by the barriers. 

4.2 Simulations 

Similar hysteresis in the motion of the drop perpendicular to the grooves is seen in lattice 
Boltzmann simulations. A liquid drop was allowed to equilibrate on the corrugated surface. 
Then, using this as an initial condition, we ran two sets of simulations where the drop volume 
was quasi-statically increased in the first set and decreased in the second set. To do this, 
we varied the liquid density by ±0.1% from its equilibrium bulk density every 10000 time 
steps. This changed the drop volume as the system relaxed back to its coexisting equilibrium 
densities. 

When the drop volume is increased, the drop contact line advances. Therefore, we expect the 
first set of simulations to mimic a drop sprayed onto the surface and the final drop shape to 
be elongated parallel to the grooves. When the drop volume is decreased, to the contrary, the 
drop contact line recedes. We therefore expect the second set to mimic a drop stamped onto 
the surface and the drop shape to be elongated perpendicular to the grooves. The simulation 
results shown in Fig. [2] (d) and (e) clearly show that the drop shape is indeed elongated 
parallel and perpendicular to the grooves when the contact line is (d) advancing and (e) 
receding. The drop contact angles were also tracked during the simulation and the results 
are shown in Fig. El The main features indicating contact line pinning in the perpendicular 
direction are again reproduced: (i) the parallel contact angle stays roughly constant and (ii) a 
large variation in the perpendicular contact angle measurements is observed. Note, however, 
that due to numerical limitations, it was only possible to follow the drop as it moved over one 
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groove so that Fig. [6] shows just one period of the repeating sawtooth pattern obtained in the 
experiments infH In the figure the receding contact line jump in the perpendicular direction 
occurs at y ~ 0.3 x 10^. 

5 Drops spreading on corrugated surfaces: understanding the 
contact angle measurements 

In this section we attempt to provide a consistent picture of the geometry of a drop spreading 
on corrugated surfaces over a wide range of surface parameters. In particular, we shed light 
on the influence of three dimensionless surface parameters: the roughness factor, the aspect 
ratio of the barriers and the relative ratio of the dimensions of the drop to the dimensions 
of the barriers. Our conclusions will be drawn from both experiments and lattice Boltzmann 
simulations. 

We have demonstrated that the shape of a drop spreading on a corrugated surface is no 
longer spherical. The drop is elongated either in the direction parallel or perpendicular to the 
grooves and the contact angle varies along the contact line. In order to simplify the analysis, 
we focus only on two principal directions: parallel and perpendicular to the grooves, since the 
greatest difference in contact angle is expected between these two. We define the apparent 
contact angle as the contact angle made by fitted ellipses at its intersection with the top 
of the ridges. It is very important to realize, however, that the apparent angle is often not 
equal to the local contact angle, in particular where the contact line is distorted locally by 
the corrugations. We also define the aspect ratio as the ratio of the groove height to the 
pitch and the drop elongation as the ratio of the maximum base radius in the parallel and 
perpendicular directions. 

5.1 Experiments 

Fig. El^a) and (b) show our experimental results on how the drop elongation and contact angles 
depend on the aspect ratio. In this set of experiments, we prepared photo-embossed surfaces 
having different corrugations and liquid drops were placed onto the surfaces with a syringe. 
We first found that there was no obvious relationship between the final drop geometry and 
either the groove depth or the groove pitch. However, plotting the data against the aspect 
ratio of the corrugations, trends could be identified, although the scatter was significant. This 
is expected given the important role of hysteresis and the existence of multiple local energy 
minima. 

There are four aspects of the data which we wish to comment on: 

(i) On chemically striped surfaces, it has been shown experimentally [3] and theoretically 
|36j that there is no contact angle hysteresis (at least due to the surface patterning) in the 
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Figure 7: Experimental results for water drops on photo-embossed surfaces, (a) Elongation of the droplet 
and (b) Static 6\\ and 9±_ as functions of the aspect ratio of the surface. Both the elongation and 6± increase 
approximately monotonically with increasing aspect ratio. The solid line corresponds to the Wenzel angle, 
(c) AS — 6± ~ S|| as a function of elongation. The solid line is the relation between A6 and the elongation 
predicted by Eq. [TOland using — ^wenzei- 
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direction parallel to the stripes. Furthermore, Gleiche et al. [1] showed that the parallel 
contact angle is close to the Cassie-Baxter contact angle |40j . Similarly, in this study, we 
might expect the parallel angle to approximately follow the Wenzel formula |41j . 



cos 9\ 



r cos 6, 



(7) 



where r is the roughness factor. In the experiments 



wi + ^Jwl +Ah? 



(8) 



V 



where ~ W2- wi, W2, h and p are respectively the ridge width, the groove width, the 
groove height, and the pitch. The equilibrium contact angle ~ 68°. In figure Wi)^) the 
Wenzel approximation is indeed shown to correspond well to ^y. The same conclusion was 
also found in [22j . 

(ii) The perpendicular contact angle, on the other hand, shows a very large scatter. Typically 
9± > and 0± tends to increase with the aspect ratio. This may be explained using the 
Gibbs' criteria [38] we discussed in the previous section. Since the drop spreads outwards 
on the corrugated surface, it is appropriate to consider the advancing contact line motion. 
Therefore, it is easier for the drop to spread parallel to the grooves. This immediately implies 
that the drop contact angle perpendicular to the grooves is larger than that parallel to the 
grooves. The maximum value of 9± is an estimate of the advancing angle in the perpendicular 
direction. For the two dimensional model illustrated in Fig. ^ 9a = 9^ + a, where a is the 
maximum slope on the surface. This is clearly an overestimate in three dimensions, since 
it neglects the energy costs associated with the surface deformation from a spherical cap. 
Nevertheless, from the data shown in Fig. \7\ we found that the advancing angle increases 
from 70° to 88° with increasing aspect ratio. This is consistent with the fact that, for the 
experimental substrate geometry, (Fig. [I] (a)), the maximum slope a increases with the aspect 
ratio of the grooves. Similar behaviour is also observed in Fig. 2(e) of a recent work by Chung 
et al. [22]. 

(iii) In Fig. [7] (c), we plot the difference between the contact angles A9 as a function of 
the drop elongation. We find that A9 varies monotonically with drop elongation and the 
experimental scatter is reduced when compared to figure [7] (a) and (b). Qualitatively this 
can be explained by assuming that the drop shape is close to elliptical. Assume for simplicity 
that the drop profiles in the two principal directions can be fitted to circles. In this case, the 
apparent contact angles can be written as 



where a and b are the drop base lengths in the two principal directions. As a result. 



tan6'||/2 = h/a 



and 



tan6'_L/2 = h/b, 



(9) 



tan6l_L/2 = e tan 6*11/2. 



(10) 



15 



50 100 150 200 250 

' 50 100 150 200 25tf 

Figure 8: A typical drop shape observed in simulations, (a) and (b) the drop cross sections in the direction 
parallel and perpendicular to the grooves, (c) Top view of the drop. The dashed lines are fits to (a) and (b) a 
circle or (c) an ellipse. 

^11 is approximately constant in the experiments and hence 6± increases monotonically with 
e. In Fig. [TJ^c), Eq. is plotted assuming = ^wenzcb and again it corresponds well to 
the trend observed in the experiments. 

(iv) The drop elongation e increases monotonically with the aspect ratio. This result follows 
from the fact that the energy gain of spreading in the direction parallel to the grooves is 
higher for higher aspect ratio. 

5.2 Simulations 

To gain further insight into these relationships between surface and drop geometry, simulations 
were performed. The drop was initialized as a spherical drop just above the surface and was 
allowed to spread, without any external force, on the corrugated surface. The drop elongation 
and contact angles were then recorded after 400000 time steps, by which time the drop had 
reached its (meta) stable configuration. A typical final drop shape is shown in Fig. [HI As 
shown in Fig. M^a) and (b) the drop cross sections in the direction parallel and perpendicular 
to the grooves can be fitted well with ellipses. We have also tried to fit the contact line around 
the base of the drop, see e.g. figure [8l^c), but the typical fit we obtain is poor: indeed, the top 
view of the drop clearly shows a corrugated contact line. The corrugation is more noticable 
with increasing aspect ratio. This is due to the fact that it is more advantageous to wet the 
sides of the grooves the higher the aspect ratio. This is similar to changing the wettability 
ratio between the hydrophilic and hydrophobic stripes on chemically striped surfaces. 

The simulation results are presented in Table [1] and Fig. [9l In interpreting the data it 
is important to note that there were two differences to the experiments forced upon us by 
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Table 1: Drop contact angles and elongation for different groove heights. The drop volume is ~ 7.47 x 10^ 
and the groove width and pitch are kept constant at 30 and 40 respectively. 6'wcnzci is the theoretical Wenzel 
angle. Here the drops only lie on 4 grooves. 

computational requirements. Firstly, the dimensions of the simulated drop were comparable 
to the dimension of the grooves so that a drop lies on only 4 to 5 barriers. Secondly, the 
ridges were taken to be rectangular so that the maximum slope a = 90° was independent of 
the aspect ratio. We varied the groove depth but kept the other parameters the same. In the 
simulations 

r = l + — . (11) 
P 

Of the experimental features listed in the previous subsection, (iii) increases with in- 
creasing elongation e, and (iv) the drop elongation increases with increasing aspect ratio, are 
reproduced in the simulations, as shown in Fig. and (c). These trends are observed for 
the same reasons as in the experiments. 

The experimental trends (i) and (ii) are, however, not reproduced by the simulations for the 
following reasons: 

(i) ~ ^Wcnzci: in the simulations (shown in Fig. [Hb)), the Wenzel angle is a poor ap- 
proximation for the parallel contact angle, even though both simulations and experiments 
are in the same aspect ratio regime. The reason for this discrepancy is that the size of the 
heterogeneities is of the same order as the size of the drop. Since the drop only lies on a small 
number of grooves, the Wenzel equation is not applicable. For this reason we performed a 
second set of simulations where the size of the ridges is an order of magnitude smaller than 
the drop, as shown in Table [5] and Fig. [TUl and found the Wenzel equation is indeed a rea- 
sonable approximation to the parallel contact angle. In this case, we note that the size of the 
ridges is comparable to the width of the interface in the lattice Boltzmann simulations and 
the roughness factor is much higher than that obtained in the experiments. In other words, 
as long as the typical length scale of the corrugation is much smaller than the drop size, the 
Wenzel equation is a good estimate of the parallel contact angle. 

(ii) 0± > 0|| and 6± on average increases with increasing aspect ratio: from Fig. [U^b) we 
indeed found that 9± > ^m, however, from the data we have so far, increasing the aspect ratio 
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Figure 9: Lattice Boltzmann simulation results for drops on a corrugated surface, (a) Elongation of the 
droplet and (b) 6\\ and 9±_ as function of the aspect ratio of the surface. The elongation increases monotonically 
with increasing aspect ratio, but both 9\\ and decrease. The straight line corresponds to the Wenzel formula, 
(c) as a function of the elongation. (The dotted lines in (a) and (c) are guides to the eye.) 
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Table 2: Drop contact angles and elongation for different groove widths. The drop volume is ~ 7.42 x 10^ 
and the barrier width and height are kept constant at 4 and 4 respectively. 6'wonzoi is the theoretical Wenzel 
angle. 
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Figure 10: and 9± as a function of the aspect ratio of the surface. Here the size of the ridges is much 
smaller than the size of the drop. The solid line is the theoretical Wenzel angle. 

is not followed by an increase in 6±; instead a decrease was found. We now argue why this 
may be the case. In the experiments, the roughness factor is close to 1 (at most 1.06, from 
equation[8]) and hence 6*11 is typically just below 6^- Furthermore, 9a increases with the aspect 
ratio since a increases with the aspect ratio. Since < 9± < 9^, 0± increases with increasing 
aspect ratio. In the simulations, on the other hand, r is considerably larger than 1 (between 
1.15 and 2) and hence is much smaller than 9e. Since the lower limit of the allowed values 
of 9± is decreased, it is plausible that we find 9± decreasing with increasing aspect ratio in 
Fig. Mh). 



6 Conclusions 

We have used experiments and simulations to investigate the behaviour of drops on surfaces 
patterned with sub-micron and micron-scale parallel grooves. We find that the final drop 
shape is highly dependent on the path by which it is achieved. Drops which advance across 
the surface are elongated parallel to the grooves whereas drops that dewet the surface are 
elongated perpendicular to the grooves. We explain this behaviour in terms of the pinning 
of the contact line on the groove edges. We stress that it is not possible to measure a single 
contact angle for a drop on a patterned surface. The contact angle varies around the rim of 
the drop, with its position on the surface and with its dynamic history. 

Contact angles and base dimensions were measured for drops spreading on corrugated surfaces. 
This was done for a large number of drops and substrate geometries. The data was very noisy, 
underlining the prevalence of hysteresis on patterned surfaces. However, we were able to draw 
several conclusions from the results. We find that the parallel contact angle is close to the 
Wenzel angle. The Wenzel angle in return depends on the roughness factor of the surface. 
Pinning effects render results for the perpendicular contact angles very noisy but they are. 
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in general, greater than the parahel angles and tend to increase with increasing aspect ratio. 
The difference between parallel and perpendicular contact angles increases with the drop 
elongation whilst the drop elongation increases with the aspect ratio. This behaviour is 
explained, again by considering contact line pinning, and by noting that the drop spreads 
more quickly in the parallel direction on surfaces with higher aspect ratio. 

We found that lattice Boltzmann simulations of drop motion were of great use in interpreting 
the experiments. The same relation between the direction of elongation of the drop and its 
direction of motion over the surface was immediately apparent in the simulations, and we 
were able to probe the pinning perpendicular to the grooves. Differences in detail between 
the simulations and experiments arose because the simulations are limited to drops lying on 
a small number of grooves. In particular, we find that the Wenzel angle no longer provides a 
good estimate of the parallel angle when the drop dimensions are comparable to the dimen- 
sions of the corrugation. An interesting future project would be to perform experiments on 
wider grooves where a quantitative match to the simulations may be possible. 
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